function result=examoupath(sigmai,rd)

a=1;m=1;dt=0.5/5000;b=1;
result=sigmai+a*(m-sigmai)*dt+b*rd*sqrt(dt);

